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ABSTRACT 

A  2D  digital  image  S  is  represented  conventionally  by  the  union  of  grid  squares 
containing  pixels  of  5  which  we  denote  by  F(S).  This  gives  the  correct  topology  for  5 
with  8-adjacency,  and  with  a  little  imagination,  4-adjacency  can  also  be  properly  han¬ 
dled.  However,  one  encounters  difficulty  in  extending  basic  2D  results  to  3D  digital  im¬ 
ages.  The  last  few  years  have  seen  the  need  for  better  methods  which  give  a  closer  link 
with  well  developed  continuous  topology,  especially  with  the  advent  of  digital  surface 
theory  [11].  We  define  a  new  continuous  model  F{S)  by  refining  F{S).  We  show  that 
this  gives  a  better  bridge  between  the  two  subjects,  digital  and  continuous  topologies. 
We  also  show  how  this  space  F{S)  is  related  to  two  other  continuous  models  [4]  [7]. 
Although  we  concentrate  only  on  2D  images  in  this  paper,  the  concepts  and  general  ideas 
extend  to  3D  images.  A  3D  version  of  this  paper  is  in  preparation. 
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1.  Introduction 


A  digital  image  is  a  discrete  sample  of  a  continuous  picture  in  real  life.  The 
sample  is  usually  taken  on  an  evenly  spaced  grid  point  set  on  a  rectangular 
domain  of  a  plane.  To  be  definite,  think  of  the  set  of  points  with  integral  coordi¬ 
nates  in  the  x-y  plane  as  the  grid  point  set. 

For  example,  suppose  the  sample  looks  like  the  figure  below: 

110  0 
0  0  11 

Figure  1.1.  The  I’s  indicate  pixels  of  an  object,  and  the  O’s  indicate  pixels 
of  the  background. 


The  natural  question  is  “Does  this  digital  sample  suggest  one  (connected) 
object  or  two  separate  ones  of  similar  appearance?”  Of  course,  it  also  depends  on 
how  the  sample  was  taken,  but  this  is  not  our  concern  here.  We  are  interested  in 
defining  geometric  concepts  on  digital  images  consistently  and  realistically,  which 
is  about  all  we  can  do  without  knowing  the  original  picture.  Returning  to  the 
connectedness  question,  if  we  regard  the  set  of  I’s  shown  above  as  a  connected 


set,  then  the  set  of  O’s  should  be  a  disconnected  set  in  order  to  be  consistent  with 
the  standard  topology  on  the  Euclidean  plane.  Thus  it  is  customary  in  digital 
topology  to  use  different  adjacency  topologies  on  a  digital  image  (denoted  by  5 ) 

•r 

and  on  its  background  (regarded  as  the  complementary  set,  and  denoted  by  5 ). 
For  example,  if  8-adjacency  is  used  on  S,  which  makes  the  I’s  in  Figure  1.1  con-  j 
nected,  then  4-adjacency  is  used  on  5  which  makes  the  O’s  disconnected. 
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In  order  to  approximate  the  unknown  original  picture  by  a  continuous  one, 
it  has  been  the  tradition  to  represent  each  point  p  in  the  digital  image  S  by  the 
closed  unit  square  around  p  ,  and  to  represent  5  by  the  union  of  all  such  squares. 
Let  us  denote  this  union  by  F  (S ).  This  gives  a  good  topological  representation 
for  the  8-adjacency  topology  on  5  in  the  sense  that  S  is  connected  in  this  topol¬ 
ogy  if  and  only  if  F(S)  is  connected  in  the  usual  Euclidean  topology.  Of  course, 
we  cannot  say  the  same  thing  about  the  4-adjacency  topology  at  the  same  time. 
However,  it  is  not  difficult  to  describe  the  connectedness  property  using  generally 
accepted  terms  of  Euclidean  topology  in  F  (S ). 


(a)  (b) 

Figure  1.2.  (a)  F(S ),  the  traditional  pixel  representation,  (b)  E  (5 ),  the  skeletal 
representation  used  for  computing  the  Euler  characteristic  number 
of  5. 

An  alternative  to  the  square  representation  has  existed  all  along,  dating  back 
to  Minsky  and  Papert’s  1969  Perceptrons  [8].  In  computing  the  Euler  charac¬ 
teristic  number  of  a  digital  Image  S ,  they  took  a  point  for  each  point  in  5 ,  a  line 
segment  for  each  pair  of  4-adjacent  points,  and  a  square  for  each  2X2  matrix  of 
four  4-adjacent  points  in  5 .  This  gives  a  good  skeletal  representation  of  S  with 
4-adjacency  which  will  be  denoted  by  E{S).  Of  course,  it  is  not  difficult  to  deal 
with  8-adjacency  as  well.  Curiously  enough,  all  subsequent  authors  (except  Mor- 
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genthaler  [10])  used  this  representation  only  in  computing  Euler  characteristic 
numbers  until  Kong  and  Roscoe  used  it  for  other  purposes  in  1985  [7|. 

As  long  as  one  is  concerned  with  only  2D  images,  there  is  no  real  difficulty  in 
defining  continuous  representations  of  digital  images.  The  situation  is  different 
when  we  consider  3D  images.  A  number  of  basic  facts  about  2D  images  were 
generalized  to  3D  images  without  great  difficulty,  leading  to  the  initial  develop¬ 
ment  of  3D  digital  topology  by  Rosenfeld  in  1980  [15],  but  one  crucial  theorem 
for  this  foundation  did  not  yield  to  an  easy  generalization.  It  is  Proposition  5  in 
[15].  The  stumbling  block  was  the  fact  that  the  corresponding  theorem  for  2D 
images  used  a  boundary  curve  tracking  algorithm  which  did  not  extend  to  surface 
tracking  in  3D  in  any  obvious  way.  Meanwhile  Herman  and  his  colleagues  were 
developing  such  a  surface  tracking  algorithm.  Based  on  topological  results  from 
the  classical  book  by  Newman  [12],  Herman  and  Webster  proved  [4]  that  their 
previously  developed  algorithm  [l]  [2]  was  correct.  However,  the  surface  traced 
by  their  algorithm  is  not  a  surface  in  the  strict  topological  sense.  Nevertheless, 
the  main  result  of  this  work  was  essential  to  the  proof  of  Proposition  5  given  by 
Morgenthaler  and  Rosenfeld  in  [ll]. 

During  the  same  period  of  time,  a  theory  of  digital  surfaces  was  developed 
[11].  A  surface  was  defined  as  a  digital  image  satisfying  a  certain  set  of  axioms 
analogous  to  the  axioms  for  a  manifold  in  topology.  Although  this  definition  is 
clever  and  simple,  it  was  no  easy  task  to  prove  the  digital  analog  of  the  standard 
fact  in  topology  that  a  closed  surface  in  Euclidean  3D  space  is  orientable  [14]. 
Thus  it  was  natural  to  seek  for  some  closer  connection  between  the  two  subjects. 
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the  newly  developed  3D  digital  topology  and  the  well  established  Euclidean  topol¬ 
ogy,  so  that  one  could  use  known  facts  from  this  latter  discipline.  Recently  Kong 
and  Roscoe  did  precisely  that  in  1985  (7).  They  defined  the  continuous  analog  of 
a  digital  image  as  the  3D  version  of  the  representation  considered  by  Minsky  and 
Papert  and  others  as  mentioned  above,  and  proved  that  a  digital  image  5  is  a 
surface  in  the  sense  of  Morgenthaler  and  Rosenfeld  [llj  if  and  only  if  its  continu¬ 
ous  analog  is  a  surface  in  the  strict  topologfcal  sense.  In  other  words,  the 
axiomatic  definition  of  a  digital  surface  given  by  [ll]  was  correct  (if  direct  adja¬ 
cency  is  iised  on  the  digital  image,  and  indirect  adjacency  is  used  on  the  back¬ 
ground  image). 

The  purpose  of  this  paper  is  to  propose  another  method  of  describing  a  digi¬ 
tal  image  5  by  a  continuous  representation  to  be  denoted  by  F  (5 ).  This 
representation  F  (S )  seems  more  consistent  with  the  topological  behavior  of 
Euclidean  space  in  the  sense  that  it  tends  to  give  easier  access  to  basic  results 
from  that  well  developed  subject.  We  first  give  a  new  proof  of  the  main  theorem 
of  Herman  and  Webster  [4]  from  a  conceptual  point  of  view.  Their  proof  is  a 
combinatorial  analysis  of  different  types  of  strong  connectedness.  We  reduce  all 
these  to  the  usual  connectedness  in  topology  using  F  (S ).  We  use  this  as  a 
motivation  for  defining  F  (S ),  but  the  advantages  of  F  {S)  should  be  apparent 
once  one  sees  the  definition.  To  define  it,  we  introduce  a  concept  of  singularity 
similar  to  that  used  in  algebraic  geometry  and  differential  topology.  We  then 
adopt  one  of  several  possible  ways  of  removing  the  singularities  from  F  {S  ).  The 
resulting  topological  space  f  (5 )  is  a  compact  manifold  with  a  boundary  in  the 


strict  topological  sense.  In  particular,  the  boundary  is  a  true  surface  in  the  3D 
case,  and  a  true  curve  in  the  2D  case.  This  surface  is  an  approximation  to  the 
surface  traced  by  the  surface  tracking  algorithm  mentioned  above. 

We  begin  with  brief  descriptions  of  F(S)  and  [/  (S ),  and  with  the  introduc¬ 
tion  of  cellular  complexes  in  a  very  restricted  sense  (Section  2).  In  the  following 
two  sections,  we  define  a  2D  analog  of  the  surface  tracking  algorithm  by  means  of 
a  graph  [1]  [2]  [4],  and  we  give  a  definition  of  F  (5 )  while  constructing  a  proof  of 
the  2D  analog  of  the  theorem  due  to  Herman  and  Webster  [4].  Combining  the 
results  in  Sections  2  and  4,  we  give  a  proof  of  Proposition  5  [ll|  (Section  5). 
Finally,  the  three  continuous  representations  of  a  digital  image  are  shown  to  be 
intimately  related  (Section  6).  For  example,  the  2D  analog  of  the  above  men¬ 
tioned  surface  tracking  algorithm  can  be  used  to  trace  the  continuous  analog 
£(S). 

The  content  of  this  paper  is  devoted  to  2D  images,  but  it  is  so  arranged  as 
to  yield  an  easy  extension  to  3D  images  in  principle. 


2.  Continuous  representations 


We  assume  that  the  reader  is  familiar  with  the  basic  definitions  and  elemen¬ 
tary  facts  about  2D  digital  geometry  given  in  [16].  Let  E  be  the  set  of  all 
integral  points  within  a  large  rectangular  domain  in  R  and  let  5  be  a  subset  of 
E.  For  each  p  in  E,  define  F  (p)  to  be  the  unit  closed  square  with  p  at  its 
center,  that  is  the  set  of  all  points  {x ,  y)  in  such  that 

I®  “  P  il  ^  ly  “  P2I  ^  Define  F{S)  to  be  the  union  of  F  (p  )  for  all  p  in  5 , 
and  U{S)  to  be  the  interior  of  F(S).  (In  analogy  to  F(p  ),  one  could  define  the 
open  set  U  (p  )  for  p  in  5  to  be  the  interior  of  F  {S ).  However,  this  would  be 
quite  useless  since  the  union  of  the  U{p)  for  p  in  5  consists  of  disjoint  open 
sets.)  Notice  that  F{S)  (resp.  C/^(5))  is  connected  in  the  sense  of  the  standard 
Euclidean  topology  if  and  only  if  5  is  8-  (resp.  4-)  connected  in  the  sense  of  digi¬ 
tal  topology.  In  other  words,  we  can  say  that: 

The  8-  (reap.  4‘)  adjacency  topology  on  a  £D  digital  image  S  ia  well  represented 
by  the  standard  Euclidean  topology  on  F{S)  (reap.  U  (5  )). 

The  topological  spaces  F{S)  and  U{S)  will  be  referred  to  as  continuous 
representations  of  S  (in  the  primitive  form).  E  and  hence  5  being  a  finite  set, 
F  {S  )  (resp.  U  {S  ))  is  a  bounded  and  closed  (resp.  open)  subset  of  R  Thus  if  5 
is  represented  by  F  {S)  then  the  complementary  set  5  (called  the  background) 
should  be  represented  by  U  {S ).  If  we  assume  that  S  does  not  meet  the  boun- 
daty  of  E,  we  have 

C/(E)  =  f(5)U  U(S). 


B 


This  is  the  basic  reason  for  using  the  8-adjacency  topology  on  S  if  the  4- 
adjacency  topology  is  used  on  5. 


Now  we  study  global  properties  of  the  topological  space  F  {S)  by  means  of 
the  adjacency  relations  among  the  constituent  blocks  F  (p )  for  p  in  5 .  The  set 
together  with  the  set  of  these  blocks  is  a  special  case  of  a  cellular  complex(*) 
which  we  shall  not  define  precisely  here.  One  relevant  remark  is  that  it  is  some¬ 
times  important  to  distinguish  the  single  set  F  {S)  from  the  collection  of  all  its 
constituent  cells:  all  squares  F  (p)  for  p  in  S ,  all  their  edges,  and  all  their  ver¬ 
tices.  If  we  denote  this  collection  by  K  then  we  denote  the  union  of  its  consti¬ 
tuent  cells  by  \K  |  which  is  F  (S ).  The  figure  below  shows  an  example, 

p  0 

(a)  9  »•  (b) 

Figure  2.1.  (a)  5  =  {p  ,q  ,r  }.  (b)  F  (S ). 

where  K  consists  of  3  squares,  10  edges,  and  8  vertices.  Fortunately  the  cell 
complex  K  made  out  of  the  cells  in  F{S)  is  uniquely  determined  by  the  single 
set  F{S).  A  cellular  l-complex(*)  is  a  collection  of  edges  in  some  F[S)  together 
with  all  their  vertices.  If  Z,  is  a  cellular  1-complex  then  \L  |  denotes  the  union  of 
all  edges  in  L  .  For  example,  let  F  =*=  F  (p  )  Pi  F  (? ),  the  common  edge  of  two 
squares,  and  let  «  ,  v  be  the  two  end  points  of  E .  Then  L  =  {E  ,u  ,v}  is  a  cel¬ 
lular  1-complex.  For  any  pixel  p  in  E,  define  the  (cellular)  boundary  dF  (p  )  to 
be  the  1-complex  consisting  of  its  four  edges  and  the  four  corner  vertex  points. 


Notice  that  the  set-theoretic  boundary  Bd  F{p)  in  R’  is  simply  |5F(p)|.  If 
E  ^  F{p)f\F{q)  with  p  in  S ,  q  in  5  then  we  denote  p  =  i{E),  q  =  o  {E) 
(i  for  inside,  o  for  outside).  Of  course,  ail  these  depend  on  5 .  We  define  the 
(cellular)  boundary  dF  {S )  to  be  the  1-complex  consisting  of  all  edges  E  and  the 
end  point  vertices  of  E  for  which  there  are  pixels  p  in  5 ,  g  in  5  with 
E  =  F  (p )  n  F  (^ ).  We  leave  to  the  reader  the  proof  of  the  following  impor¬ 
tant  but  not  difficult  lemma: 

Lemma  \dF  (5  )\  =  Bd  F  (5 ). 


(*)The  cellular  complex,  also  known  as  the  CW  complex,  was  initially 
defined  by  the  English  mathematician  J.H.C.  Whitehead  in  the  1930’s.  It  con¬ 
sists  of  very  general  types  of  cells.  It  was  invented  for  studying  the  deeper  struc¬ 
ture  of  a  topological  space,  such  as  multiple-connectedness.  It  is  used  for  com¬ 
puting  global  properties.  For  example,  the  Euler  characteristic  number  of  a  CIV 
complex  can  be  computed  by  counting  the  number  of  cells  in  each  dimension 
which  is  shown  to  depend  only  on  the  underlying  topological  space  independent 
of  any  particular  cellular  complex  on  the  space.  See  Spanier  (1966)  [17]  for  a 
more  comprehensive  exposition  on  this  subject.  Of  course,  we  need  only  consider 
the  very  special  case  arising  from  digital  topology.  Klette  formulated  such  a  spe¬ 
cial  type  of  cellular  m -complex  in  Euclidean  n -space.  (See  Klette  (1983)  [6].)  It 
would  be  useful  to  establish,  in  this  setting  suitable  to  digital  topology,  elementa¬ 
ry  facts  analogous  to  those  in  algebraic  topology.  The  only  work  involving 
homology  groups  in  digital  topology  appears  in  Janos  and  Rosenfeld  (1981)  [5]. 
This  work  too  touches  very  lightly  upon  an  algebraic  method  from  topology. 
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3.  Graph  representation 


In  this  section  we  shall  define  a  2D  version  of  the  graph  due  to  Herman  et  al. 
[l]  [2]  [4]  as  a  data  structure  for  a  boundary  curve.  Following  their  notation,  let 
Q  —  F{S)  he  a,  cellular  complex,  and  let  P  =  R“  -  Q .  Two  edges  E ,  F  in 
dQ  are  said  to  be  P-adjacent  if  {!)  E ,  F  are  adjacent,  i.e.,  E ,  F  have  a  com¬ 
mon  vertex,  and  (2)  o  {E )  =  o  {F ),  or  o{E),  o{F)  are  4-adjacent,  or  o  {E ), 
0  (F )  are  4-adjacent  to  a  common  pixel  in  S.  Notice  that  the  first  condition 
determines  a  2X2  matrix  containing  all  of  the  pixels  o{E),  i{E),  o  {F ),  i  (F ), 
and  that  the  second  condition  is  equivalent  to  the  statement  that  (2^)  o  {E ), 
o{F)  are  4-connected  within  the  2X2  matrix.  There  are  six  possible  cases  of 
relative  positions  of  o  {E ),  i  {E ),  o  {F  ),  and  i  {F )  within  the  2X2  matrix: 


E  E 
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1 

1 

0 

0 

1 

1 

1 
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0 

1 

0 
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0 

1 

0 
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0 

0 

1 

0 

0 

(a)  (b) 


(c)  (d) 

Figure  3.1. 


(e)  (f) 


The  edges  E,  F  in  Figures  3.1  (a),  (b),  (c),  (f)  are  P-adjacent  whereas  those  in 
the  other  two  figures  are  not.  A  P -component  of  dQ  is  defined  as  usual.  Define 
a  graph  for  dQ  by:  For  each  edge  in  dQ .  create  a  node  n  {E ),  and  for  each  P- 
adjacent  pair  of  edges  E ,  F  ,  connect  the  nodes  n{E),  n{F  )  by  an  arc.  We  shall 
denote  the  graph  by  'y{Q  ).  We  can  turn  )  into  a  digraph  by:  For  each  p  in 
E,  orient  F  {p  )  clockwise,  and  let  the  induced  orientation  on  the  four  edges  be  as 


indicated  below; 


< - 

Figure  3.2. 


Thus  if  =  F  (p  )  Pi  F  (7  )  then  the  two  orientations  on  E  induced  by  those  on 
F  (p  ),  F  {q)  are  opposite.  For  each  edge  E  in  9Q ,  we  shall  use  the  orientation 
on  E  induced  by  0  {E ).  Thus  E  is  a  directed  line  segment,  written  E  —  ui. 
Define  an  ordered  pair  E ,  F  to  be  P-adjacent  if  the  end  point  of  the  directed 
line  segment  of  E  is  the  beginning  point  of  the  directed  line  segment  F ,  that  is, 
if  F  =  ut)  then  F  =  Ctfi.  Clearly  P-adjacency  implies  F -adjacency.  A  quick 
survey  of  all  six  cases  in  Figure  3.1  shows  that  if  F ,  F  are  P-adjacent  then 
either  F ,  F ,  or  F ,  F  are  P-adjacent.  This  then  makes  it  possible  to  direct 
every  arc  of  the  graph  n(Q  }•  We  shall  denote  the  resulting  digraph  by  r(Q). 
The  examples  below  reveal  some  interesting  facts  about  these  graphs. 


1] 


Figure  3.4.  (a)  5.  (b)  Q  =  F(5).  (c)  dQ  with  directed  edges,  (d)  r(Q  ). 


Notice  that  in  Figure  3,  r(Q)  has  three  components  whereas  dQ  has  only  two, 
and  that  in  Figure  4,  both  dQ  and  r(Q )  have  only  one  component  each. 
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4.  Continuous  representations  in  the  refined  form 

The  purpose  of  this  section  is  to  introduce  a  concept  of  singularity  in  a  digi¬ 
tal  inaage  S,  and  to  show  how  to  remove  the  corresponding  singular  point  from 
the  topological  space  F{S)  to  construct  a  better  continuous  model  for  S .  We 
shall  do  this  while  deducing  a  2D  version  of  the  main  theorem  in  Herman  et  al. 
[4]  from  a  purely  topological  theorem  proved  in  Appendix  A. 

Let  Q  F  {S)  and  P  ~  -  Q  as  before.  Let  C  be  a  component  of  Q , 

and  let  Z)  be  a  component  of  R'^  ~  C . 

Proposition  4.1.  Bd  D  is  a  P-component  of  dQ  . 

In  general,  a  P-component  is  not  a  component  of  dQ  as  seen  from  the 
example  in  Figure  3.3c,  and  a  P  -component  is  not  a  simple  closed  curve  as  seen 
from  the  example  in  Figure  3.4c.  It  is  clear  from  the  definition  that  P- 
components  are  connected  sets. 

We  need  only  prove  that  (1)  every  pair  of  points  of  Bd  D  are  in  the  same 
P-component;  (2)  if  F  C  dQ  is  an  edge  in  the  same  P-component  of  an  edge  E 
in  Bd  D  then  F  C.  Bd  D .  First  notice  that  Bd  D  is  a  cellular  1-complex  in 
which  each  edge  E  has  the  property  that  U {o  {E  ))  G  D  ,  F  {i  {E  ))  C.  C .  For  (2) 
we  need  only  prove  it  when  F  is  actually  P -adjacent  to  E,  but  then  o{E), 
0  (F  )  are  in  the  same  4-component  in  P  ^  -  C  and  » {E  ),  i  [F  )  are  8-connected  in 
Q ,  and  hence  F{i{F))  C  C,  U{o{F))G  D  implying  F  G  Bd  D  as  claimed  in 


For  (1),  it  follows  from  Proposition  A  in  Appendix  A  that  Bd  D  is  con¬ 
nected.  Thus  if  we  know  that  components  and  P -components  are  the  same  then 
we  are  done.  A  study  of  ail  six  cases  of  adjacent  pairs  of  edges  E ,  F  \n  dQ  in 
Figure  2  showed  that  non-P -adjacency  occurred  only  in  a  2X2  square  of  the 
0  1 

form  j  Q.  Such  a  square  (or  its  rotation  by  90  degrees)  shall  be  referred  to  as  a 

singularity  in  S .  The  non-existence  of  a  singularity  in  S  guarantees  that  any 
adjacent  edges  in  dQ  are  P -adjacent.  Thus  if  we  assume  that  there  is  no  singu¬ 
larity  in  S  then  we  are  done. 

0  p 

Now  assume  the  existence  of  singularities  in  5 .  If  ^  0  ^  singularity  then 

the  common  vertex  point  P(p)n  P(?)  will  be  called  a  singular  point  in  F{S). 
The  key  problem  in  representing  S  by  P  (5 )  arises  from  the  presence  of  singular 
points:  F  {S )  looks  like  ^  ^  near  a  singular  point  x  and  is  connected  via 
a  single  point.  At  the  same  time,  the  background  set  near  x  consists  of  two  dis¬ 
joint  open  sets  but  they  are  so  close  that  their  boundaries  meet  though  only  at  a 
single  point.  This  problem  can  be  solved  by  modifying  F  (5 )  slightly  near  each 
singularity  of  5 .  Incidentally,  this  kind  of  problem  exists  in  many  other  fields 
such  as  algebraic  geometry,  differential  topology,  etc.  There  are  several  possible 
ways  of  resolving  singularities,  but  here  we  shall  adopt  the  simplest  one  which 
may  be  achieved  in  practice  by  taking  a  sufficiently  find  grid.  The  basic  idea  is 
to  widen  the  bottle-neck  at  x  from  F  (p )  to  F  (g  )  without  altering  the  global  pic¬ 
ture  of  F(S).  For  each  singular  point  z,  simply  add  to  P  (5 )  a  l/4-size  closed 
square  with  its  center  at  z . 


o 


Figure  4.1.  (a)  Near  a  singular  point  x  in  F  {S ).  (b)  The  singular  point 
X  resolved  in  F  (5  )  by  adding  a  small  square  around  x  . 

Let  F  (5 )  =  F(5)  (^J,-  cr^  )  where  i,  ranges  over  all  singular  points  of 

F  (5 ).  One  could  use  a  unit  size  square  for  ,  but  this  would  alter  the  shape  of 
F(S)  too  much.  Now  F(5)  has  no  singular  points.  Write  Q  =  Q  U  (U> 

C  =  C  \J^  (ijy  )  with  ly  in  C ,  ^  =  D  -  [Jy  .  Then  C  is  a  component 
in  Q  and  D  is  a  component  in  In  this  case,  Bd  Dt  has  just  been 

proved  to  be  a  component  oi  dQ  .  Hence  Bd  Dt  is  a  simple  closed  polygon,  and 
the  corresponding  Bd  D  is  a,  F -component  of  dQ  .  // 


The  space  F  {S )  gives  more  accurate  topological  descriptions  of  the  digital 
image  S ,  and  its  border  sets  with  respect  to  8-adjacency  on  5 .  The  alteration  of 
F  (5 )  altered  U{S)  as  well.  Define  C/  (5 )  to  be  U  (E)  -  F  (5 ).  We  shall  refer  to 
F(5),  U{S)  as  continuous  representations  of  S  with  8-adjacency  fin  the  refined 
form  ). 

We  can  define  F  (5 )  analogously  for  4-adjacency  on  S ;  in  this  case  the  open 
square  around  x  is  removed  from  F  {S  )  for  each  singular  point  x  in  F  (5  ).  Thus 
F  (5 )  depends  on  the  type  of  adjacency  used  on  S .  When  there  is  possibility  of 
confusion,  we  shall  use  F^  (5 )  for  x -adjacency. 


TO 


Morgenthaler  considered  the  problem  of  finding  conditions  on  S  under  which 
the  topological  space  F(S)  correctly  represents  the  digital  topology  on  S.  The 
conditions  he  proposed  were  phrased  in  terms  of  local  connectedness.  It  turns  out 
that  those  conditions  were  devised  to  guarantee  the  non-existence  of  a  singularity 
in  S .  We  shall  prove  this  in  Appendix  B. 

Before  closing  this  section,  we  wish  to  pont  out  that  resolution  of  singularity 
can  be  achieved  on  the  digital  level.  First  subdivide  the  grid  space  into  half-units 
on  each  axis.  Thus  each  unit  square  is  subdivided  into  4  smaller  squares. 
Replace  each  pixel  by  4  new  pixels  of  the  same  kind.  So  far,  there  is  no  real 
change;  a  singularity  remains  a  singularity.  Now  we  replace  each  singularity  in 
5  as  in  the  figure  below: 
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Figure  4.2.  (a)  A  singularity,  (b)  Finer  grid,  (c)  Singularity  resolved  for 
8-adjacency  on  S .  (d)  Same  for  4-adjacency  on  S . 


It  may  be  wiser  to  use  a  finer  grid  for  a  closer  geometric  approximation.  We 


actually  used  a  1/4  unit  length  size  grid  in  the  definition  of  F  (S ). 


5.  Application  to  a  connectedness  theorem  in  digital  topology 

The  following  theorem  on  2D  images  was  proved  in  Rosenfeld  and  Kak  [16] 
and  its  3D  analog  was  proved  indirectly  using  the  surface  tracking  algorithm  due 
to  Herman  et  al.  in  [1]  [2]  [4].  The  purpose  of  this  section  is  to  show  how  the 
method  of  continuous  representation  of  the  previous  section  may  be  used  in 
deducing  this  basic  theorem  in  digital  topology  from  the  topological  result  of  the 
previous  section. 

Proposition  5.1.  Let  S  be  a  digital  image.  Let  C  be  a  4-component  of  S, 
and  D  an  8-component  of  S  such  that  C,  D  are  4-odjacent.  Then  Cp  is  4- 
connected  where  Cp  is  the  set  of  p  in  C  which  are  8-adjacent  to  D. 

Let  us  first  prove  the  following  exact  digital  analog  of  the  theorem  of  the 
previous  section: 

Lemma  5.1.  Let  S  be  a  digital  image.  Let  C  be  an  8-component  of  S,  and 
let  D  be  a  4-component  of  C  .  Then  Dq  is  4-connected  where  Dq  is  the  set  of  all 
p  in  D  which  are  8-adjacent  to  C. 

Let  Q  =  F{S).  Then  F{C)  is  a  component  of  F(S).  Assume  that  S  does 
not  meet  the  border  of  S  and  that  U  (E)  has  no  boundary  so  that  F  {S)  G  U  (E). 
Take  any  pair  p,  q  in  Dc-  We  want  to  show  that  p,  q  are  in  the  same  4- 
component.  By  assumption,  p  is  8-adjacent  to  C .  Assume  p  is  not  4-adjacent 
to  C .  Then  there  is  a  pixel  a  in  C  located  relative  to  p  forming  a  2X2  matrix 

p  X  _ 

^  Here  the  other  two  pixels  ’s)  must  be  in  C.  They  are  4-adjacent  to  p. 
and  hence  are  in  the  same  4-component  as  p  ,  namely  D .  Let  b  be  one  of  these 


X ’s.  Then  the  edge  E  —  F  (a  )  Pi  F (6  )  is  in  fid  D  .  Similarly  there  are  pixels 
c  in  C ,  d  in  D  such  that  c  ,  d  are  4-adjacent  and  q ,  d  are  also  4-adjacent. 
The  edge  G  =  F  (c  )  Pi  F  (d )  is  in  fid  D  .  Now  by  Proposition  4.1,  E ,  G  are 
in  the  same  F -component.  This  implies  that  b ,  d  are  in  the  same  4-component 
in  Dq  ,  and  hence  p  ,  q  are  also.  / / 

To  prove  the  proposition,  let  D'  be  the  unique  8-component  of  C  containing 
D .  In  view  of  the  above  lemma,  we  need  only  show  Cp  =  Cpi.  The  relation 
“C”  follows  immediately.  Observe  that  Cp  is  the  set  of  pixels  p  in  C  such  that 
F  (p )  Pi  5F  (fi )  0.  Similarly  for  Cpi.  This  reduces  to  proving 
dF  [D')  C  dF  {D  ).  Remember  that  dF  {D  )  is  the  union  of  all  edges  E  such  that 
E  =  F  (p  )  Pi  F  («? )  for  p  in  D ,  q  in  D.  Similarly  dF{D')  is  the  union  of  all 
edges  E'  such  that  E'  =  F(p')  Pi  F(9')  for  p'  in  D',  q'  in  D'.  Actually  q’  Is 
in  C  in  the  above  expression  of  E'.  It  q'  ^  C  then  q'  €  D'  since  p',  q'  are  4- 
adjacent,  and  D'  is  an  8-component  of  C .  We  need  only  prove  that  if  E ,  E’  are 
adjacent  edges  such  that  E  C  dF{D)  and  E'  C  dF[D')  then  E'  C  dF[D). 
Write  E'  =  F(p')  PI  F(9')  as  explained  above.  Now  we  need  only  prove  that 
p'  E.  D  and  q'  E  D  .  Since  C  Pi  fi  =0,  q'  S.  C  implies  q'  E  D  .  To  prove 
q'  E  D,  consider  all  possible  cases  of  E,  E':  If  p  =  p'  then  p'  E  D  since 
p  E  D  .  Suppose  X  or  y  is  p  in  the  figure  below; 


X 

y 

P’ 

q’ 

E’ 
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If  p'  E  S  then  p'  €  C  since  p\  q’  are  4-adjacent,  which  is  absurd.  Thus  we 
must  have  p'  E  S,  but  then  p'  E  D  since  p'  is  8-adjacent  to  p.  Finally  notice 
that  p  =  q'  \s  not  possible  since  p  E  D  C.  S  and  q'  E  C  C.  S .  j j 
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6.  Relations  among  various  continuous  representations 

In  this  section,  we  shall  describe  how  closely  the  three  different  continuous 
representations  of  a  digital  surface  are  related.  For  a  2D  digital  image  5  with  4- 
adjacency,  define  a  cellular  complex  E{S)  whose  vertices  are  the  mid-points  v  (p  ) 
of  the  squares  F  (p  )  for  p  in  5 ,  whose  edges  are  the  line  segments  connecting  4- 
adjacent  pairs  in  S ,  and  whose  squares  are  the  convex  hulls  of  the  four  vertices 
of  2X2  matrices  of  4-adjacent  pixels.  Notice  that  this  is  slightly  more  general 
than  the  definition  given  in  Section  3.  The  example  below  should  clear  up  any 
ambiguity  about  this  definition. 

0  0  10  0 
110  11 
0  0  111 

(a)  (b) 

Figure  6.1.  (a)  5.  (b)^(5). 

One  can  easily  extend  this  definition  to  higher  dimensions  with  various  types  of 
adjacencies.  Kong  and  Roscoe  did  an  extensive  study  of  3D  digital  topology  in 
[7],  in  which  they  prove  that  a  digital  image  5  with  6-adjacency  is  a  digital  sur¬ 
face  in  the  sense  of  Morgenthaler  and  Rosenfeld  [ll]  if  and  only  if  the  3D  analog 
E  (5 )  is  a  surface  in  the  strict  topological  sense.  We  shall  prove  in  essence  a  2D 
analog  of  the  theorem:  If  a  3D  digital  image  5  is  a  digital  surface  in  the  sense  of 
[ll],  then  F{S)  is  topologically  equivalent  to  the  cross  product  E{S)  X  [0,1] 
{E{S)  fattened),  to  be  denoted  by  E(S )  where  [0,1]  is  the  closed  unit  interval. 


Moreover,  the  surface  tracking  algorithm  can  be  used  to  trace  any  boundary  com¬ 
ponent  of  dF{S),  and  hence  any  component  of  E(5 ).  For  any  pixel  p,  denote 
by  N{p)  the  set  of  all  8-neighbors  excluding  p .  Recall  that  a  pixel  p  in  S  is  a 
4-curve  point  if  the  following  2D  analog  of  a  simple  surface  point  is  satisfied: 

(p  )  n  5  has  exactly  two  components  4-adjacent  to  p  ,  iV  (p  )  Pi  S  has  exactly 
two  components  4-adjacent  to  p,  and  every  q  \n  N{p)C\S  is  8-adjacent  to 
both  of  these  components.  5  is  a  digital  4-curve  if  every  pixel  of  S  is  a  4-curve 
point.  This  defines  a  digital  analog  of  a  compact  1-manifold  without  boundary  in 
the  sense  of  continuous  topology.  Such  a  space  is  known  to  be  a  finite  family  of 
simple  closed  curves.  A  compact  connected  1-manifold  with  non-empty  boundary 
is  known  to  be  an  arc.  Therefore  a  compact  1-manifold  is  a  finite  union  of  a  dis¬ 
joint  family  of  arcs  and  simple  closed  curves.  For  a  2D  digital  image,  the  follow¬ 
ing  equivalent  definition  will  be  more  convenient:  A  simple  closed  curve  is  a 
sequence  of  pixels  Pi,  P2>  ■  ■  ■  ,  Pn  >4)  such  that  p,-  is  4-adjacent  to  Pj  if  and 
only  if  »  =  y  -f-  1  (mod  n ).  The  restriction  n  >4  is  to  exclude  the  two  degen¬ 
erate  cases:  S  consisting  of  a  single  pixel,  and  S  consisting  of  a  single  2X2 
1  1 

matrix  ^  i*  digital  4-arc  is  defined  similarly. 

The  following  lemma  is  a  2D  version  of  the  analogous  theorem  due  to  Kong 
and  Roscoe  [7]. 

Lemma  6.1.  A  digital  image  S  is  a  simple  closed  4-curve  (resp.  4-urc)  if 
and  only  if  E{S)  is  a  simple  closed  polygon  (resp.  arc). 


Proof.  Let  5  be  a  simple  closed  4-curve.  Using  the  second  definition 
above,  it  is  clear  that  E{S)  does  not  contain  an  isolated  point  or  a  square,  and 


that  if  S  consists  of  p^,  p  2,  ■  ■  ,  p„  then  E{S)  consists  of  vertices  v(p,  )  and 

line  segments  joining  only  every  consecutive  pair  v  (p, ),  v  {pj )  for 
j  =1-1-1  (mod  n  ).  In  other  words,  E{S)  is  a  simple  closed  polygon. 

Conversely  let  E{S)  be  a  simple  closed  polygon  with  vertices 
Vi,  V  2,  ■■■  ,  Vn  labeled  in  consecutive  order.  Let  p,  be  the  unique  pixel  in  S 
such  that  V  (p, )  =  v,- .  Since  an  isolated  single  point  is  not  regarded  as  a  simple 
closed  curve,  S  cannot  contain  just  one  pixel.  Similarly  S  cannot  be  a  single 
1  1 

2X2  matrix  like 

Hence  we  must  have  n  >4.  Now  it  is  clear  that  the  sequence 
P  v  P2>  ■  ■  ■  >  Pn  satisfies  the  condition  for  a  simple  closed  curve.  (The  proof  for 
the  statement  in  the  parentheses  is  similar.)  / / 

The  following  theorem  gives  a  partial  characterization  of  a  simple  closed  4- 
curve,  and  at  the  same  time,  it  reveals  a  close  relation  between  the  two  continu¬ 
ous  representations  of  a  digital  image. 

Proposition  6.2.  Let  S  be  a  2D  digital  image.  If  S  is  a  simple  closed  4- 
curve  then  F  {S)  —  E(5  )  (homeomorphism). 

Proof.  Suppose  that  S  is  a  simple  closed  4-curve.  If  S  contains  no  singu¬ 
larity  then  F  {S)  ^  E  {S),  and  F  {S )  is  clearly  homeomorphic  with  E(5 )  as 
claimed,  and  we  are  done.  We  study  how  a  singularity  may  occur  in  S .  Pick 
any  pixel  p  in  5.  By  definition,  there  are  exactly  two  other  pixels  4-adjacent  to 
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p .  There  are  only  two  possible  cases  (up  to  an  equivalent  trcinsformation)  as  in 
the  figure  below: 


XXX 

XXX 

1  p  1 

1  p  X 

XXX 

X  1  X 

(a) 

(b) 

Figure  6.2. 


No  singularity  can  occur  in  Figure  6.2(a).  A  singularity  may  occur  in  Figure 
6.2(b),  but  only  in  the  upper  right  2X2  sub-matrix.  This  shows  that  F{S)  is 
definitely  not  a  product  with  [0,1]  but  that  F  {S)  is  as  claimed.  / / 

The  converse  is  not  true.  Let  5  be  a  digital  set  of  two  4-ad.’acent  pixels. 
Then  F (5 )  as  F (5 ),  and  certainly  F{S)^  E{S),  but  5  is  not  a  simple  closed 
4-curve.  Notice  that  5  is  an  arc.  Indeed  we  have 

Proposition  6.3.  //  F  (5 )  ss  E(5 )  (homeomorphism)  then  S  is  a  disjoint 
union  of  ./-arcs  and  simple  closed  ^-curves,  and  conversely. 

Proof.  Assume  F  (5 )  Si  E(5 ).  There  are  one-one  correspondences  between 
the  following  sets:  components  of  E(5  ),  components  of  E  (5 ),  4-components  of 
S ,  components  of  F  (5 ).  The  conclusion  follows  from  Lemma  6.1. 

The  converse  follows  easily.  / / 

So  far,  a  digital  analog  of  a  surface  with  boundary  has  not  been  defined  in 
the  literature.  With  a  suitable  definition,  one  should  be  able  to  prove  the  follow¬ 
ing  conjecture:  A  3D  image  S  is  a  surface  with  boundary  if  and  only  if 


F  {S)  ^  E(5 )  (homeomorphism). 

Remark  6.4.  The  2D  analog  of  the  surface  tracking  algorithm  of  Herman 
et  al.  traces  dF  {S ). 

A  rigorous  proof  is  left  to  the  interested  reader.  We  shall  merely  study  the 
generic  examples  given  in  Figures  3.3  and  3.4.  Let  S  be  the  digital  image  with 
4-adjacency  in  Figure  3.3a  which  is  a  simple  closed  4-curve.  The  F  (S )  is  a  fat¬ 
tened  simple  closed  polygon  E  {S ).  The  suggested  tracking  algorithm  can  be 
used  for  S  with  8-adjacency,  and  it  traces  out  the  boundary  curves  of  4- 
components  of  the  complement  of  5 ,  that  is  just  S .  In  our  example,  there  are 
two  separate  boundary  curves.  The  interior  curve  looks  like  Figure  3.4c,  which  is 
topologically  equivalent  to  its  graph  in  Figure  3.4d.  It  is  an  important  fact  that 
the  curves  traced  by  the  above  mentioned  tracking  algorithm  are  not  in  general 
simple  closed  polygons. 

Next  let  5  be  the  same  image,  but  with  8-adjacency.  Then  there  are  three 
components  in  U  (S ),  two  of  which  are  holes  and  the  third  of  which  is  outside  5 . 
The  previously  stated  algorithm  can  be  applied  to  Q  =  F  (S ).  It  traces  out 
three  P -components  of  dQ .  In  our  example,  they  look  like  Figure  3.3c  with  the 
two  interior  squares  separated  from  each  other.  Topologically,  they  look  more 
like  the  graph  in  Figure  3.3d  which  is  just  dF(S).  (Caution:  The  space  F  (S ) 
depends  on  the  adjacency  used  on  5.)  // 

In  this  sense,  we  can  say  that  the  algorithm  of  Herman  et  al.  traces  E  {S ). 


7.  Concluding  remarks 


We  defined  the  closed  subspace  F  (S )  and  the  open  subspace  U  {S)  of 
Euclidean  space  as  primitive  approximation  of  a  2D  digital  image  S  by  con¬ 
tinuous  topological  spaces.  This  space  F  (S )  is  not  new.  In  fact,  it  has  always 
been  thought  of  as  the  set  of  pixels  representing  5 .  We  pointed  out  that  F  {S ) 
(resp.  U{S))  is  a  correct  model  for  S  itself  with  8-  (resp.  4-)  adjacency,  but  not 
for  its  borders.  However,  one  can  get  around  this  difficulty  if  one  works  with  2D 
Images  only,  but  not  if  we  want  to  develop  3D  digital  geometry  as  witnessed  by 
the  difficulty  encountered  in  the  proof  of  orientability  of  a  (closed)  digital  surface 
(with  6-adjacency),  and  by  the  fact  that  no  direct  proof  of  Proposition  5  in  [lO] 
has  appeared  in  the  literature  up  to  now  (♦).  The  existing  proof  [11]  depends  on 
the  continuous  model  of  surfaces  defined  by  Herman  and  his  colleagues  [1]  [2]  [4]. 
Thus  it  is  clear  that  a  good  continuous  model  is  needed  so  as  to  enable  us  to 
make  use  of  the  wealth  of  results  from  continuous  topology.  In  this  paper,  we 
refined  the  model  F  (S )  to  give  a  more  accurate  description  of  S  and  its  borders. 
We  first  defined  the  concept  of  singular  points  in  F  (5 ),  and  showed  a  procedure 
to  resolve  them  to  obtain  the  refined  model  denoted  by  F  (5 ).  We  produced  a 
new  proof  of  the  main  theorem  of  Herman  and  Webster  [4]  using  F  (5 )  by  first 


(*)  We  believe  that  the  reason  for  this  is  that  digital  images  are  defined  on  a 
simply  connected  space.  All  existing  proofs  of  Proposition  5  depend  on  the 
theorem  of  Herman  and  Webster  whose  proof  is  based  on  topological  results  on 
polyhedral  sets  in  Euclidean  space,  which  is  more  than  simply  connected.  In  a 
private  conversation  with  T.Y.  Kong,  we  agreed  that  no  digital  analog  of  Euclide¬ 
an  space  has  ever  been  effectively  used  so  far.  Therefore,  if  anybody  wishes  to 
pursue  digital  geometry  from  a  purely  logical  standpoint,  it  will  be  necessary  to 
further  develop  the  digital  analog  of  Euclidean  geometry. 


constructing  a  purely  topological  analog,  and  then  examining  the  differences 
between  the  digital  and  the  analogous  topological  statements.  We  gave  another 
proof  of  Proposition  5,  but  still  depending  on  the  result  of  Herman  and  Webster 
[4]. 

Traditionally,  there  has  existed  another  continuous  model  of  S  which  we 
denoted  by  E  (S ).  It  is  roughly  U  (5 )  thinned  down.  Strangely,  E  {S )  has 
always  been  used  only  for  computing  the  Euler  characteristic  number  of  S ,  until 
the  appearance  of  Kong  and  Roscoe  (7)  in  1985.  (See  Minsky  and  Papert,  Per- 
ceptrons  (1969)  [8],  and  Rosenfeid  and  Kak  (1980)  [16].)  In  [7],  it  is  proved  that  a 
3D  digital  image  5  is  a  surface  with  6-adjacency  in  the  sense  of  Morgenthaler 
and  Rosenfeid  [11]  if  and  only  if  the  3D  analog  of  E{S)  is  a  surface  (without 
boundary)  in  the  sense  of  continuous  topology. 

We  showed  that  if  S  is  2D  analog  of  a  digital  surface  in  the  sense  of  [11] 
then  F  (5 )  is  topologically  equivalent  to  E(5 )  (the  cross  product  of  E{S)  and 
the  unit  interval).  The  converse  is  in  general  not  true,  but  if  we  modify  the 
statement  by  replacing  “surface”  by  “surface  with  boundary”  then  it  may  be 
true.  As  a  corollary,  the  2D  analog  of  the  surface  tracking  algorithm  [1]  [2]  may 
be  used  to  trace  E  (S )  in  principle. 


References 

1.  E.  Artzy  and  G.T.  Herman,  Boundary  detection  In  3-dimensions  with  a  med¬ 
ical  application,  Computer  Vision  Graphics  Image  Processing  15,  92-123 
(1981). 

2.  E.  Artzy,  G.  Frieder,  and  G.T.  Herman,  The  theory,  design,  implementation 
and  evaluation  of  a  three-dimensional  surface  detection  algorithm.  Computer 
Vision  Graphics  Image  Processing  15,  1-24  (1981). 

3.  C.R.  Dyer,  Computing  the  Euler  number  of  an  image  from  its  quadtree. 
Computer  Vision  Graphics  Image  Processing  13,  270-276  (1980). 

4.  G.T.  Herman  and  D.  Webster,  A  topological  proof  of  a  surface  tracking  algo¬ 
rithm,  Computer  Vision  Graphics  Image  Processing  23,  162-177  (1983). 

5.  L.  Janos  and  A.  Rosenfeld,  Digital  connectedness:  An  algebraic  approach, 
Pattern  Recognition  Letters  1,  135-139  (1983). 

6.  R.  Klette,  m-dimensional  cellular  spaces,  TR-1281,  Computer  Vision  Labora¬ 
tory,  Computer  Science  Center,  University  of  Maryland,  College  Park,  MD, 
May  1983. 

7.  T.Y.  Kong  and  A.W.  Roscoe,  Continuous  analogs  of  axiomatized  digital  sur¬ 
faces,  Computer  Vision  Graphics  Image  Processing  20,  60-86  (1985). 

8.  M.  Minsky  and  S.  Papert,  Perceptrons:  An  Introduction  to  Computational 
Geometry,  MIT  Press,  Cambridge,  MA  (1969). 

9.  D.G.  Morgenthaler,  Three-dimensional  digital  topology:  The  genus,  TR-980, 
Computer  Vision  Laboratory,  Computer  Science  Center,  University  of 

2S 


Maryland,  College  Park,  MD,  November  1980. 

10.  D.G.  Morgenthaler,  Three-dimensional  simple  points:  Serial  erosion,  parallel 
thinning,  and  skeletonization,  TR-1005,  Computer  Vision  Laboratory,  Com¬ 
puter  Science  Center,  University  of  Maryland,  College  Park,  MD,  February 
1981. 

11.  D.G.  Morgenthaler  and  A.  Rosenfeld,  Surfaces  in  three-dimensional  digital 
images.  Information  and  Control  51,  227-247  (1981). 

12.  M.HA.  Newmann,  Elements  of  the  Topology  of  Plane  Sets  of  Points,  second 
edition,  Cambridge  University  Press,  1961. 

13.  C.M.  Park  and  A.  Rosenfeld,  Connectivity  and  genus  in  three  dimensions, 
TR-156,  Computer  Vision  Laboratory,  Computer  Science  Center,  University 
of  Maryland,  College  Park,  MD,  May  1971. 

14.  G.M.  Reed  and  A.  Rosenfeld,  Recognition  of  surfaces  in  three-dimensional 
images,  Information  and  Control  53,  108-120  (1982). 

15.  A.  Rosenfeld,  Three-dimensional  digital  topology.  Information  and  Control 
50,  119-127  (1981). 

16.  A.  Rosenfeld  and  A.C.  Kak,  Digital  Picture  Processing,  2nd  edition. 
Academic  Press,  New  York  (1982). 

17.  E.H.  Spanier,  Algebraic  Topology,  McGraw-Hill,  New  York  (1966). 


APPENDIX  A 


The  following  is  a  topological  analog  of  the  main  theorem  of  Herman  and 
Webster  (4]. 

Proposition  A.  Let  K  be  a  connected  cellular  complex  in  ,  and  let  D 
be  a  component  of  R"  -  K  .  Then  the  boundary  set  Bd  {D  )  is  connected. 

(A  “cellular  complex”  here  can  be  thought  of  as  a  finite  union  of  cubical 
cells  as  defined  in  Section  2,  or  more  generally  as  a  CVK -complex  as  defined  in 
Spanier  [I7j.) 

We  need  a  formula  from  algebraic  topology  involving  the  Betti  numbers. 
For  a  topological  space  K,  a  sequence  of  numbers  called  Betti  numbers  6,  (T)  are 
defined  for  »  =  0,  1,  •  •  •  .  60(1')  is  the  number  of  (path)  components  of  Y,  and 
6  i(K)  is  the  number  of  independent  loops  in  Y.  For  example,  b  ^(Y)  =  2  (or  Y 
a  tonis  with  one  hole,  and  6  i(  F )  =  0  for  Y  =  /?  "  .  We  need  the  following  for¬ 
mula:  60  (-^  ^i(-^  +  where  A  ,  B  are 

both  open  subspaces  of  i?  "  with  A  B  or  A  ,  B  are  both  subcomplexes 
of  a  cellular  complex  with  A  C\  B  ^  This  formula  can  be  derived  from  the 
Mayer- Vietoris  exact  homology  sequence  for  the  pair  {A  ,  B).  (See  Spanier  [17] 
for  details.) 

Proof.  It  follows  easily  that  R'*  -  K  has  finitely  many  components,  one  of 
which  is  D .  Denote  the  other  components  by  D,-  (i  =  1,  2,  ■  ■  ■  ).  Then 
i?"  -K  -D  u(u,  A  ).  Notice  that  each  D,  is  open  in  R  "  since  X  -  K  is 
open  and  iZ"  is  locally  connected.  Let  K*  =K  U  ((J,  A)-  Then 


D  =  R”  -K*,  and  K*  —  R*  -  D  Is  a  closed  set  of  i?  " .  D  being  open  in 
/? "  ,  Z)  n  Bd  Z)  =  0.  Hence  Bd  D  C  Bd  K .  Similarly  Bd  Df  C  Bd  K  for 


all  t . 

Claim:  K*  is  connected.  It  suffices  to  prove  that  if  is  a  non-empty  sub¬ 
set  of  Zii*  which  is  open  and  closed  in  K*  then  A  =  K* .  Notice  that  A  must 

be  closed  in  R "  since  K*  is.  K  being  connected,  either  A  Pi  ZiC  =  0  or 

K  C.  A  .  Case  1:  Suppose  A  Pi  /(T  =0.  Then  A  Pi  Z),-  0  for  some  i .  Z), 

being  connected,  Z),-  C  A .  On  the  other  hand,  A  being  closed  in  Z?  ” , 
C1(Z),- )  C  A  and  hence  Bd  Z),-  =  C1(Z>, )  -  Z),-  C  A  ,  but  Bd  Z),  G  K ,  contrad¬ 
icting  the  assumption  A  Pi  Zif  =0.  Case  2:  Suppose  K  C  A  .  If  Z),  Pi  A  ^9 
then  Di  G  A  as  observed  before.  Hence  A  =  ZC  Dj )  where  j  ranges 

over  a  subset  of  the  t ’s.  Suppose  Z>,-  Pi  A  =  0.  Let  x€  Bd  Z),  G  K .  Since 
C1(Z?,' )  C  ZiC*  and  K  G  A  and  A  is  open  in  K* ,  A  is  an  open  neighborhood  of 
X  in  K* .  This  then  means  A  Pi  Z?,-  ^0  contradicting  the  assumption 

Z?,-  Pi  A  =  0.  Hence  A  must  contain  all  Z?,  ’s,  which  means  A  =  K* .  This 

proves  the  claim. 

Now  we  apply  the  formula  for  A  =  C1(Z? ),  B  =K*.  Notice 
Bd  D  =  A  r\B  andZ2’‘=AUB,  Since  C1{Z?  )  =  boiK*)  =  1,  and 
^0  (Z2  "  )  =  0,  we  conclude  that  bQ{Bd  D)  =  1.  j  j 

This  proposition  is  not  in  the  most  general  form,  which  is  of  no  interest  to 
us.  The  important  question  here  is  “What  makes  the  topological  form  so  easy 
whereas  the  purely  digital  form  seems  so  difficult?”  The  critical  difference  seems 
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to  be  the  lack  of  efficient  use  of  the  fact  that  digital  images  are  taken  in 
Euclidean  space.  The  above  proof  uses  the  property  of  simple  connectedness  of 
Euclidean  space.  Perhaps  that  is  ail  we  need  from  Euclidean  space.  If  we  replace 
by  a  non-simply  connected  space,  the  proposition  may  not  be  true.  We  end 
this  section  with  a  brief  explanation  of  this  phenomenon. 

A  space  Y  is  said  to  be  simply  connected  if  Y  is  (path-wise)  connected,  and 
every  loop  in  Y  is  deformable  to  a  point.  For  example,  an  open  n-ball  and  an 
n -sphere  in  are  simply  connected,  but  a  simple  closed  curve  is  not.  The 

above  theorem  is  not  true  if  we  replace  “simply  connected”  by  the  weaker  condi¬ 
tion  “connected”.  For  example,  take  an  annulus  for  X  as  in  the  figure  below, 
and  K  for  the  right  half  shown  by  the  shaded  area.  Then  D  is  the  left  half  in 
the  figure,  and  the  set  Bd  D  consists  of  two  disjoint  (vertical)  line  segments 
which  contradicts  the  conclusion  of  the  theorem.  This  is  because  X  is  not  simply 
connected. 
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APPENDIX  B 


Morgenthaler  postulated  conditions  on  S  under  which  S  is  a  correct  sample 
of  F  (5 )  [9].  It  turns  out  that  his  conditions  are  equivalent  to  non-existence  of 
singularities  in  S .  Denote  by  N‘{p  )  the  set  of  all  8-nelghbors  of  p  including  p , 
by  O  (5 )  the  number  of  components  in  the  sense  of  S ,  and  by  O  {X)  the  number 
of  components  of  a  topological  space  X .  His  conditions  are: 

Ml.  0{F{S  r\N{p)]  =  0(5  nyV(p)),  all  p  in  5. 

M2.  0(F(5  niV(p))  =  0(5  nAr(p)),ailp  in  5. 

This  gives  us  an  opportunity  to  characterize  singularity  from  a  different  per¬ 
spective. 

Lemma.  The  following  conditions  on  a  2D  digital  image  5  are  equivalent  to 
each  other: 

51.  5  has  no  singularity. 

52.  If  p  ,  q  in  S  are  8-adjacent  to  each  other  then  there  is  r  in  S  which  is 
4-adjacent  to  both  p  and  q  . 

55.  Analog  of  S2  for  5  . 

1  0 

Proof.  If  S2  is  false  then  there  exists  a  2X2  matrix  q  ^  (or  its  rotation  by 
90  degrees)  which  is  a  singularity.  This  proves  that  Si  implies  S2.  By  symmetry, 

51  implies  S3.  If  Si  is  false  then  5  contains  a  singularity  which  is  simply  a  2X2 
matrix  such  as  that  above  which  certainly  violates  both  S2  and  S3.  Thus  each  of 

52  and  S3  implies  Si.  // 
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Proposition.  Condition  (Ml  and  M2)  is  equivalent  to  non-existence  of  a 
singularity  in  S . 

Proof.  Recall  that  O  ( /’  [  5  O  N\p  ) )  |  equals  the  number  of  8- 

components  of  S  D  N'(p ).  Thus  if  8-adjacency  is  used  on  S  then  Ml  holds 
automatically,  and  M2  is  easily  seen  to  be  equivalent  to  SI  as  4-adjacency  is  sup¬ 
posedly  used  on  5.  By  symmetry,  if  4-adjacency  is  used  on  5,  then  Ml  is 
equivalent  to  S2.  // 
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